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THE USE OF LIDAR FOR STRATOSPHERIC MEASUREMENTS* 

M. Patrick McCormick 
Langley Reseenrch Center 


ABSTRACT 


This paper reviews stratospheric measurenrents possible with ground-based, 
8Lirbome, and satellite-bome lidar systems. The instrvusents , basic equations, 
and formats normally used for various scattering and absorption phenomena 
measurements are presented including a discussion of elastic, resonance, 

Raman, and fluorescence scattering techniques. 


INTRODUCTION 


The heart of a lidar system is the laser which was first conceived by 
Professors Schalow and Townes in December 1958 at Columbia. It was not until 
May i960 that Dr. Theodore Maiman at Hughes Research Laboratory first reduced 
the laser principle to practice. Actually the term lidar, an acronymn for 
light detection and ranging, came before the advent of the laser. Middleton 
and Spilhaus in their 1953 book on "ifeteorological Instruments" first suggested 
the term "lidar" for light detection and ranging when discussing ceilometry. 

The honor of the first paper in which data were described using a lidar goes 
to Fiocco and Smullin for their September 1963 paper "Detection of Scattering 
Layers in the Upper Atmosphere (60-lit0 km) by Optical Radar." This was 
followed by a paper by Fiocco and Grams (May 196 k) entitled "Observations of 
the Aerosol Layer at 20 km by Optical Radar." During this early period, 
papers were being presented and published in which the potential uses of 
lasers and lidar in meteorology were being expounded. As an example, Goyer 
and Watson at NCAR published a paper on "The Laser and Its Application to 
Meteorology" in September 1963* %ron Ligda gave a paper in September 196U 
at the World Conference on Radio Meteorology in which he pointed out that 
others were also discussing the possible uses of lasers in meteorology. These 
included Dick Schotland, Dave Atlas, P. S. Carter, and R. B. Battelle. Also 
described in this presentation was the first Stanford Research Institute (SRI) 
lidar (Mark I) which was built in January 1963. Initial atmospheric returns 
were obtained in July 1963, and they were the first to apply the lidar 
technique to meteorology in the lower troposphere. 

Most lidar measurements have utilized elastic scattering techniques but more 
recently have also included inelastic scattering and absorption studies. 
Presently, many groups are working with Ildar in an attempt to characterize 
the atmosphere. For an idea of the scope of this work, see the various 


^Invited Paper: 8th International Laser Radar Conference, Drexel University, 

Philadelphia, Pa., Jtone 6-9, 19T7. 


Conference Abstracts of the Conferences on Laser Radar Studies of the 
Atmosphere. A comprehensive review of lidar measurements of the atmosphere 
can be found in E. D. Binkley (1976). 

In general, laser radars operate in the following manner. A Q-switched laser 
emits a pulse of nearly monochromatic light approximately 30 nsec in duration 
into the atmosphere. Molecules and suspended particulate matter (aerosols) 
scatter and/or absorb this radiation as the p\ilse propagates through the 
atmosphere. A small portion of this light is scattered directly back toward 
the laser. A receiver composed of mirrors and/or lenses collects this back* 
scattered radiation and directs it onto a photodetector irfiose output is 
measured as a function of elapsed time after laser emission or, therefore, 
range. The backscattered energy incident on the photodetector is examined 
spectrally at or near the laser output wavelength with color filters, inter- 
ference filters, or spectrometers. This enhances the signal -to-noise ratio 
by reducing unwanted background radiation and determines whether elastic or 
inelastic techniques are utilized. 

An exanq>le of a system for stratospheric measurements is the Langley Research 
Center’s U8" laser radar system. It consists of two temperature-controlled 
lasers (ruby and neodymium-doped glass) mounted alongside an f/lO Cassegrainian- 
configiured telescope consisting of a 48-inch diameter f/2 all metal primary 
and a 10-inch diameter secondary. The output from the detector package is 
recorded by a high-speed data acquisition system. Analog signals are 
amplified and bandwidth limited, digitized at a 5- or 10-megahertz rate with 
8-bit accuracy, and recorded on magnetic tape. Pulse count data is amplified, 
discriminated, counted at a 200 megahertz rate, and also stored on magnetic 
tape. Altitude resolution is obtained by using the variable 1-, 5-, or 
10-microsecond bin widths that are available. A 16K word storage computer is 
used to control the data acquisition system and provide data massage. An 
X-band microwave radar, boresighted with the laser system axis, is used to 
ensure safe operation in the atmosphere. A rotating shutter in front of the 
laser reauces laser fluorescence after Q-switching. The entire system is 
mobile and can scan in elevation and azimuth at a slew rate of 1° per second. 


LIDAR THEORY 


After crossover of the emitted laser pulse and receiver field of view, the 
voltage V (X,R) at the photodetector output for a laser radar system is 
given by 


V(X,R) 


y(X )E(X)q^(X,AR)q^U' ,AR)f(X,R) 


U) 


with 


2 


f 


y(X) » “ T(X)S(X)L 


where f(X,R) is the backscattering function at wavelength X and range R, 
and ®^e the transmissivities from the laser radar to 

the scattering volume and from the scattering volume to the laser radar, 
respectively. The symbol X' refers to a wavelength that may not be X. L 
is the photodetector load resistor, T(X) is the system optical efficiency, 

S(X) is the photodetector spectral sensitivity, A is the area of the receiver, 
c is the speed of light, and E(X) is the laser output energy per pulse. 

For elastic scattering, the emitted and return wavelengths are Identical. 

Equation (l) reduces to 


v(A R) = 

R^ 


where f(X,R) = fM(X,R) + fA^(X,R). The subscripts M and A refer to 
molecular and aerosol, respectively. The Raman scattered return is shifted 
in wavelength from the laser output and equation (l) becomes 


y(X’ )E(X) q^(X,AR)qg(X* ,AR)f (X' ,R) 


where fj.(X',R) is the Raman scattering function at R and is independent of 
aerosol scattering. The shift in wave numbers is characteristic of the 
particular molecular species in the scattering volume. Raman scattering can 
occur for €uiy incident wavelength but its scattering function is approximately 
three orders of magnitude less than molecular scattering functions. 


For the case of fluorescence scattering, X' may or may not be X with the 
scattering function being that of fluorescence scattering. It can only occur, 
however, when X corresponds to an absorption line or band in the atmosphere 
and, like Raman scattering, the return radiation at X* is characteristic of 
the fluorescing molecules. Although the cross sections for fluorescence are 
very large compared to Raman scattering, there are Inherent problems that must 
be overcome in order to utilize this technique. It now appears that the differ 
ential absorption laser radar technique (DIAIi) might offer the best hope for 
measuring pollutants with laser radar. 


An atmospheric scattering model for the elastic backscattering function at two 
laser wavelengths is shown in figure 1. The details of its derivation are 



given ‘y McCormick (l9Tl). The molecular contribution is calculated by using 
the familiar Rayleigh scattering eq.uation and the U. S. Standard Atmosphere 
(1962). The aerosol contribution is calculated assuming: Mie theory, a 

Junge r"^ size distribution from 0.125 pm to 10 pm; a l.S refractive index; 
and Rosen's (1966) particle sampling data from 5 to 26 km with a sea level 
value of 450 particles/om^ and an exponential interpolaiioti between sea level 
and 5 km. This model includes this interpolation since the aerosol munber 
density in the first 5 km is highly variable. From equation (2) it is evident 
that the laser radar signsJ. is dependent upon a combination of both aerosol 
and molecular scattering through the transmissivity, 

q^(X,R) = exp (-2 t) = exp [-2/^ 6(X,R)dR] (U) 

where t is the optical depth and R is the extinction coefficient, and 
through the aerosol and molecular scattering functions. The molecular 
scattering contribution is very easy to calculate (Rayleigh scattering) but 
the aerosol contribution cannot be easily (and not uniquely) determined. 
UsTially spherical aerosols of an assumed size distribution and refractive 
index are iised with Mie scattering theory to calculate the aerosol 
contribution. 

A combination of Raman and elastic scattering, however, might also be used 
to determine eierosol scattering functions. The ratio of the Raman return 
at two altitudes and Z2 can be written from equation (3) as 


Z^V(X' ,Z^) q(X,0-Z^)q(X' ,Z^-0)o^N(Z^) 

Z^VlX' ,Zg) q(X,0-Z2)q(X' ,Z^-0)a^n{Z^) 


where is the Raman cross section for the particular gas molecule r.f 
number density N and ci(X,0-Zx) is the transmissivity at X from the 
laser radar to altitude Z]_. 

Assuming q{X,AZ) = q(X’,AZ), equation (5) becomes 


Z^V(X’,Z^) 1 .l(Z^) 

Z^VCX'.Z^) q^(X',Z^-Z2) ^^^2^ 


( 6 ) 


This assumption is good when there is no differential absorption or large 
extinction between X and X’. Taking the natural logarithm of both sides 
of equation (6) gives 




/ 2 $(X' , Z)dZ - -fin -5 == 

Z^ Z^V(X',Zj)N(Z2) 


Assuming 8 constant over Z^ to Z^ giver 


2(Z2-Z^) 


z:v(x* Z„)N(ZJ 


Z^V(X',Z^)H(Z2) 


The number densities for air at the two altitudes can be obtained from 
rawinsonde temperatxire and pressure data. Then 3 can be used in 
Koschmleder’s equation to find meteorological range. Or this value of the 
total extinction coefficient can be used with equation {k) to determine the 
two-way transmissivity and the value of q^ obtained in this manner can be 
used in equation (2) for elastic scattering. The Raman nitrogen data ratioed 
at three altitudes combined with simultaneous elastic scattering data ratioed 
at the same three altitudes will yield three equations and three unknowns 
from which the aerosol scattering function at the three altitudes can be 
uniquely determined. Other techniques for using laser radar data are discussed 
in a paper by McCormick and Fuller (1973). 

Normally, stratospheric aerosol lidar data is presented in terms of the 
backscattering ratio 


R'(X,Z) = Z^V(X,Z)/q^fj^(X,Z) 


or using equation (2) 


R’ (X,Z) * yE{X) 


f^(X,Z) + 

f„(X,Z) 


The lidar total backscatter ratio R(X,Z) can be determined by multiplying 
R'(X,Z) with a conversion constant K such that 


R(X,Z) = K R'(X,Z) 


= 1 + f.(X,Z)/f (X,Z) 


Since the Ildar system calibration factor is unknown, the constant K has 
to be determined by normalizing the ratio in an altitude region where 

backscatterlng is assumed to be only from molecules. This normalization can be 
performed by searching for the minimum value of F*(1,Z) in the altitude range 
from the tropopause to 32 km for each profile. The exact altitude where 
aerosol backscatterlng is assun»d negligible depends on varied factors such 
as the approximate aerosol content of the stratosphere (fresh volcanic, back- 
ground volcanic, etc.), seasonal variations in the aerosol profile, lidar 
wavelengths, etc. The minimum value of R(A,Z) found is assumed to be due to 
molecular backscatterlng alone and used to normalize the backscattered profile. 
As described above the two-way atmospheric transmittance must also be determined 
either by assused theoretical models or by a combination of Raman and elastic 
lidcnr backscatterlng. The sensitivity of q2(x,Z) to various stratospheric 
aerosol loadings has been c’iscussed in detail by Russell et al. (1976). If 
the proper normalization and transmittance are assumed, any ratio vsilue 
greater than one indicates aerosol scattering. 


STRATOSPHERIC MEASUREMENTS 


Ground-Based 

Lidars have been used extensively to study the lower stratosphere. Some of the 
earlier measm*ements are shown in figure 2 taken over Williamsburg, Virginia, 
March l8, 196? (McCormick, 196?). At this period of time the stratospheric 
aerosol layer (Jxuige or sulfate layer), centered at approximately l8 km, was 
very evident in lidar backscatter signatures representing approximately twice 
the scattering expected from a molecular atmosphere. 

Shown in figure 3 are examples of this scattering ratio plotted as a function 
of altitude for the period 1967-197^. After 1967, the stratospheric aerosol 
loading decreased ar*d became quiescent until the eruption of the volcano 
de Fuego in Guatemala (l^*.5° N, 90.9® N). Large Increases in stratospheric 
aerosol and their related optical effects were observed at different locations 
in the northern hemisphere (McCormick and Fuller, 1975; Fegley and Ellis, 

1975, Melnel and Meinel, 1975; Volz, 1975). Additional data covering the 
period 196it throiagh 197^ axe shown in figure h (Russell et al., 1975). 

They are compared with three other techniques used during this period. Despite 
the fact that the four observational methods observe different characteristics 
of the stratospheric aerosol, and the locations are different, the general 
trends all agree. 

The data from another comparison study corroborating the lidar technique are 
shown in figure 5 (Northern et al., 197^). In this study the normalized lidar 
backscattering was used to infer the aerosol scattering function. The lidar 
data were compared with the Wyoming dust sonde derived aerosol concentration. 

The best fit was for aerosol scattering functions of 6 to 8 x 10~9 m~L sr~^. 

The agreement between the results from these two instruments is generally 
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quite good. Both measurements detected thin high altitude cirrus clouds, both 
describe the peaX and approximate structure of the 20 km aerosol layer and 
both indicate a relatively aerosol free region above the troposphere. 

Another comparison is shown in figure 6 (Russell et al., I97!>). Data from 
the SRI lidar is compared with filter sample data taken aboard a WB-5T aircraft 
that flew over the lidar site. A calculation based upon the mass of collected 
particles and assumed particle model parameters for a modified gamma distribu- 
tion of sulfuric acid particles yields the agreement in aerosol backscattering 
coefficient indicated. 

A more recent comparison is shown in table 1 (Pepin et al., 1976). Here, lidar 
data are compared with dustsonde measurements and satellite limb extinction 
measurements. The results from these three different measiirement techniques 
agree very well. 

Son« lidsur data obtained after the de Fuego eruption are shown in figure 7 
corresponding to four different months; November 197*»» February 1975, July 
1975, and J\xly 1976. The dashed lines are rawinsonde-derived temperature 
profiles and the arrow indicates the location of the assumed tropopause. A 
summary of the gradual deceiy of the strc*ocpneric aerosol after the de Fuego 
eruption is further shown in figure 8 (McCormick et al., 1977), where the 
aerosol backscattering coefficient profile has been integrated with altitude 
over different altitude intervals. 

28 

/ (R(X,Z)-1) • fj^(X,Z)dZ (12) 


where hi is the height of the tropopause in km. This integral has been 
normalized to the value obtained from the lidar data of October 10, 197^. 

These data show an e- folding decay time of 11.6 months over the time period 
Jsuiuarj' 1975 to July 1976. From both winter to summer periods of 1975 and 
1976 the data show an approximate 1*0 percent decrease in the aerosol integrated 
backscatter. The decrease of the total stratospheric aerosol integrated 
backscatter over a 1-year period from January 1975 to January 1976 is approximately 
65 percent. These values are in relative agreement with estimates of the 
annual outflow of stratospheric air mass from the northern hemisphere 
(Reiter, 1975 ). 

Another technique for measuring stratospheric aerosols is to use two 
wavelengths simultaneously, one sensitive primarily to aerosol scattering 
and fr.e second to molectilar scattering. Efforts have been made to analyze 
lidar backscattering data independently of rawinsonde data since the real 
time molecular scattering is needed in order to properly interpret the aerosol 
scattering ratio and the integrated aerosol column data, especially in times 
cf low stratospheric loadings. Otherwise, rawinsonde date must be collected, 
used to construct the expected molecular scattering profile, and the lidar 
ratio fcrmed as in equation (9). Figure 9 (Fuller et al. , 1976) shows lidar 
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data at 0.69^3 ym ratloed to 0.3^72 ym Ildar data as compared to Ildar data at 
0.69^3 ym ratloed with the standard ravlnsonde technique. Another approach 
to this problem is the use of Raman backscattering from nitrogen. Kent et al. 
(1971) in Jamaica measured Raman backscattering from nitrogen to a Uo km 
height. A comparison with rawinsonde-derlveu molecular density profiles showed 
good agreement. 

Inabu (1976) has suialytically shown the feasibility for measuring various gases 
in the atmosphere using Raman backscatter. His results indicate that Raman 
backscatter at night using a laser energy of 1 Joule per pulse at 0.3^72 ym 
and a receiving telescope of 0.^ m^ area,, can be used to probe oxygen emd 
nitrogen to approximately 50 km b'^lght with 100 laser shots. His simxilatlons 
use 1 km altitude intervals. H2O, N2O, NO, and O3 were below the detection 
limits in the stratosphere, but he concluded that CO2 should be measurable to 
approximately 15 km. 

Othax’ techniques such as near-resonance Raman scattering and resonance 
fluorescence are not feasible for stratospheric species measurements using 
present state-of-the-art instrumentation. Stratospheric ozone measurements, 
however, are on the threshold of present detection limits using dlfferenticd 
wavelength and range absorption lldeur measurement techniques. A balloon- 
borne or airborne lidar system for measuring various stratospheric species 
at short ranges, however, appears feasible. 


Aircraft Measurements 

Another stratospheric application of lidar was shown in a paper by Stanford 
(1975) on the temporal variation of lidar scattering from stratospheric 
condensation particles over polar regions. He proposed that these measure- 
ments could shed valuable light on the physics and meteorology involved in 
their formation by yielding information on the height, time, and horizontal 
variation of the water vapor mixing ratio. 

The ability of an airborne lidar system for probing the stratosphere was 
clearly shown in a paper by Fernald and Schuster (1977). It showed that 
aerosols could vary systematically with the prevailing meteorology and that 
continuous airborne lidar measurements over large geographical areas could 
provide valuable means for studying stratospheric /tropospheric exchemge 
processes. 


Satellite Measurements 

A number of lidar groups including SRI, Wisconsin, and NASA Langley Research 
Center are evaluating the feasibility of measuring atmospheric species from 
satellite altitudes (Wright et al. , 1975; Shipley et al., 1975; and Remsberg 
and Northern, 1976). These studies show the possibility of measuring a number 
of atmospheric species from satellite altitudes. Because of the power and 
instrument size and weight requirements, these measurements will have to be 
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perfomwd from space shuttle type platforms. The measurements that are 
feasible from such a system include: Spatial variations of cloud tops; 

subvisible cirrus cloud distribution, and in soim instances cloud thickness; 
stratospheric aerosol distributions; tropospheric aerosol distributions vdien 
not obscured by thick clouds; atomic vapors such as sodium; several minor ions 
such as Mg'^; and possibly molecular constitutents such as O2, N2, H2O, and O3. 

Figure 10 shows a stratospheric simulation using various aerosol models from 
Shettle and Fenn (1976): Fresh volcanic, aged volcanic, and background 

condition. The calculations are for a backscatter return at a spacecraft 
altitude of 260 km. The ordinate shows the munber of photons per microsecond 
Incident on a 1-meter-dlameter mirror for a shuttle-borne Nd:YAG or Nd:glass 
laser emitting a single 1 Joule pulse at I.06 ym. It is possible with this 
signed level to detect the stratospheric aerosol to at least 30 km in night- 
time viewing conditions. Figure 11 shows an exas^le of the return from a 
fresh volcanic aerosol with and %rithout a 1 km uniform cloud having an 
extinction coefficient of 2.9 km"^. Cloud top determinations are possible 
to a resolution of the Ildar r^nge cell length in all viewing conditions 
(daytime and nighttime). 
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TABLE I AEROSOL «>TICAL NOKL 



RATIO OF LUAR BACKSCATTER FIRKTION AT 0.7 TO EXTIRaiOR AT 0.83 |M. 
^OUSTSOHDE RESULTS 


Scattering function, f ( \.i ), m stersd 



rigure !•- Backscattering function model of the atmosphere for a wavelength 
of 0*69^3 pm and 0 » 3 ^T 2 pm. Molecular contribution, aerosol 
contribution, and their sum are shown. 
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(d) AIRCRAFT SAMPLE 
HEKLA 


FERNANDINA®^ 


FUEGO 


1862 1964 1966 1868 1970 1972 1974 

U.- Comparison of stratospheric aerosol measurements made by four independent techniques 
between 1962 and 197^. (a) Maximum value of peak aerosol backscatter ratio, (R-1.0) 

as observed by GF(Grams and Flocco), CL(Collis and Ligda), C(Clemesha et al . ) » 
S(Schuster), CR(Clemesha and Rodriques), O(Ottway), FS(Fresh and Schuster), F(Fox 
et al.), YE(Young and Elford), RVHC(Russell , Veezee, Hake and Collis). (b) Number 
of particles above tropopause as measured by photoelectric particle counter 
(Hofmann et al.). (c) Color-ratio an 5 )litude as measured by twilight photometer 

(Volz), (d) Mass of sulfate per standard cubic meter (scan), as collected on 
aircraft-borne filters (Castlemsui) . 













BALLOON -SONDE AEROSOL COUNTER 

OPTICAL RADAR AEROSOL MEASUREMENTS 
USING REAL ATMOSPHERIC MOLECULAR 
NUMBER DENSITY 

SEPT. 13-14 1972 


SMOOTHED AEROSOL SCATTERING FUNCTION, f, 


AEROSOL CONCENTRATION (NO./cm’) 


Figure 5*- Comparison of dustsonde particle concentration and lidar scattering 
function (9/13-1^/72) • 











Ki^'ure Y.- Lidar backscatter profiles: Post-Fuego. The dashed lines are 

rawinsonde derived temperature profiles and the aiTov indicate 
the location of the assumed tropopause. 



stratospheric layers from October 197U to July 1976 








Fifmre 9.- Lidar backscatter ratios 
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